One-step multi-qubit GHZ state generation in a circuit QED system 
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We propose a one-step scheme to generate GHZ states for superconducting flux qubits or charge 
qubits in a circuit QED setup. The GHZ state can be produced within the coherence time of the 
multi-qubit system. Our scheme is independent of the initial state of the transmission line resonator 
and works in the presence of higher harmonic modes. Our analysis also shows that the scheme is 
robust to various operation errors and environmental noise. 



Entanglement is the most important resource for quan- 
tum information processing. Therefore, the question of 
how to prepare maximally entangled states, i.e., the GHZ 
state, or the Bell states in the two-qubit case, in vari- 
ous systems remains an important issue. Superconduct- 
ing Josephson junction qubits are one of the promis- 
ing solid-state candidates for a physical realization of 
the building blocks of a quantum information proces- 
sor, see e.gji«2»2»£, They are undergoing rapid develop- 
ment experimentally, in particular, in circuit QED se- 
tups. Two-qubit Bell states have been demonstrated ex- 
perimentally^^. There are also some theoretical pro- 
posals on how to generate maximally entangled states 
for two or three qubits^ ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 . However, how 
to scale up to multi-qubit GHZ state generation remains 
an open question. Some general schemes based on fully 
connected qubit network is proposed but no specific cir- 
cuit design is provided 1 ^. Most recently, preparation of 
multi-qubit GHZ states was proposed based on measure- 
ment^'^. This type of state preparation is probabilis- 
tic and the probability to achieve a GHZ state decreases 
exponentially with the number of qubits. In this pa- 
per, we propose a GHZ state preparation scheme based 
on the non-perturbative dynamic evolution of the qubit- 
resonator system. The preparation time is short and the 
preparation is robust to environmental decoherence and 
operation errors. 



I. THE COUPLED CIRCUIT QED SYSTEM 

The GHZ state preparation scheme described below 
is based on a circuit QED setup where superconduct- 
ing qubits are strongly coupled to a ID superconducting 
transmission line resonator (TLR). FigureQJa) shows the 
type of circuit we have in mind: a qubit array is placed 
in parallel with a line of length Lq. The superconducting 
transmission line is essentially an LC resonator with dis- 
tributed inductance and capacitanc e) 20 ' 21 . The oscillating 
supercurrent vanishes at the end of the transmission line 
and this provides the boundary condition for the electro- 
magnetic field of this on-chip resonator. The qubits are 
fabricated around the central positions x = Lq/2. Since 
the qubit dimension (several micrometer) is much smaller 
than the wave length of the fundamental electromagnetic 
modes (centimeter), the coupling between the qubits and 
the TLR is approximately homogeneous. Since x — Lq/2 
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FIG. 1: (Color online) Schematic diagram of our setup, (a) 
The qubits are coupled through a superconducting stripline 
resonator (the blue stripe). Each 'crossed box' denotes one 
qubit which can be either a charge qubit or a flux qubit; the 
dashed red line shows the magnitude of the magnetic field, 
(b) Detailed schematic of a charge qubit. The crosses denote 
Josephson junctions; (c) Detailed schematic of a gradiometer- 
type flux qubit. The crosses denote Josephson junctions. 



is an antinode of the magnetic field where the electric 
field is zero, the qubits are only coupled to the magnetic 
component, which induces a magnetic flux $' through 
the superconducting loop given by 



$' = r ^)££( a + a t) 

7T 



with 



n 



(i) 



MWtt 



(1) 



(2) 



Here, AfW is the mutual inductance between the res- 
onator and the i-th qubit, u = ^/(LC) 1 / 2 is the fre- 
quency of the fundamental resonator mode, <&o = h/2e 
is the magnetic flux quantum and L (C) is the total self- 
inductance (capacitance) of the stripline. Here, we have 
assumed the qubit array to be only coupled with a sin- 
gle mode of the resonator, and a (aJ) is the annihilation 
(creation) operator of this fundamental mode. 

The stripline resonator can be used to couple both 
charge qubits and flux qubits as described below. 
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A. Charge qubit system 



lations are 



We first consider the charge qubit case. Suppose each 
qubit is a charge qubit (see Fig. QJb)) consisting of a dc- 
SQUID formed by a superconducting island connected to 
two Josephson junctions. The Coulomb energy of each 
qubit is modified by an external bias voltage and the ef- 
fective Josephson tunneling energy is determined by the 
magnetic flux $^ threading the dc-SQUID. The Hamil- 
tonian of a single charge qubit reads^ 2 - 



E. 



(i) 



3> 



(0 



H d) = ±SL (1 _ 2n W )(T W _ E f cos(7r^)a« , (3) 



where E^ (Ej) is the Coulomb (Josephson) energy of 

the i-th qubit, and rig is the bias charge number that 
can be controlled by an external gate voltage. The Pauli 
matrices a z = |0) (0| - |1) (1| , a x = |0> (1| |1) (0| are 
defined in terms of the charge eigenstates |0) and |1). |0) 
and |1) denote and 1 excess Cooper pair on the island 



respectively. $ 



$'W includes contributions 



from both the external flux bias and the flux <&'W . 

For small r/W , the Josephson energy can be expanded 
to linear order in rf 1 ^ , which results in an additional linear 
coupling between the x-component of the qubits and the 
bosonic mode. If all the qubits are assumed to be biased 
at the degeneracy point n g — 1/2, the total Hamiltonian 
reads 



H = 



(fiW($«)4 i )+ 5 «($W)(a + at) ( 7« 



H 



LC 



with the single charge qubit energy splitting fiW(<&^) = 

- J B ( ; ) cos(7r$i l) /$o), 9 {l) (^e y ) = r7 W ^ ;) sin(7r$i 4) /$o), 
and the free Hamiltonian of the TLR Hlc = ua^a. Note 
that the coupling between the qubits and the TLR can 
be turned off by setting <&e = n$o. 



B. Flux qubit system 



For a flux qubit system, a circuit example to realize our 
proposal is shown in Fig. Die). The i-th qubit contains 
four Josephson junctions in three loops instead of one or 
two loops in the conventional flux qubit desig n 23 i 24 . The 
two junctions in the dc-SQUID have identical Joseph- 
son energies affi Ej\ here is the ratio between the 
Josephson energy of the smaller junction and that of the 
two bigger junction o 23 ' 24 . The other two junctions are 

(i) 

assumed to have the Josephson energy Ej . The super- 
conducting loops are penetrated by magnetic fluxes $^ , 

(i) (i) 

! an d &d respectively. The corresponding phase re- 



(5) 



^+^+^+4- = a *(*g-*$)/*o (6) 



n$ , 



(7) 



where <pk (k = 1, 2, 3, 4) is the phase difference across the 
k-th junction. The total Josephson energy of the circuit 
is 

-U^ = Ef caaip^+Ef cos iff 

+ a^Ef cos (27r$f V*o " {<P? + (4°)) (8) 

with $f } = - and «W = 2a£ ] cos(7r$^ ) /$o)- 

If is biased close to $o/2, the circuit becomes a flux 
qubit, i.e., a two-level system in the quantum regim e 23 ' . 
Together with the charging energy, the total Hamiltonian 
for the i-th qubit is 



H [i) 



5 W(*W )£r « 



a«>(*<; 



(9) 



The Pauli matrices read a z — |0)(0| — |1)(1|, a x = 
|0)(1||1)(0|, and are defined in terms of the classical 
current where |0) and |1) denote the states with clock- 
wise and counterclockwise currents in the loop. The 
energy spacing of the two current states is eW($j ) = 
Ip\^[^ — $o/2), and the tunneling matrix element be- 
tween the two states is AW ') = AW(qW). Note that 
in contrast to the original flux qubit desig n 23 ' 24 , this gra- 
diometer flux qubit is insensitive to homogeneous fluctu- 
ations of the magnetic flux 2 ^ 5 . More importantly, it en- 
ables the TLR to couple with the dc-SQUID loop without 
changing the total bias flux of the qubit. As in the case 
of the charge qubit, the magnetic flux in the dc-SQUID 
loop includes two parts: <5>^ = ^ + where 
is due to the external control line and is due to the 
TLR. 

For »7« < 1, one can expand the Hamiltonian in terms 
of 77. The second-order terms ~ r]^ 2 d 2 A/du 2 are much 
smaller than the zeroth and the first-order term. The 
Hamiltonian of each qubit can be written as 2 ^ 

ffW = e «($« ) CT W + A«($W) ( r( i ) +5 (i)($«) CT W( a+a t) . 

(10) 

The coupling coefficient is 



0«($W) = -2af^ l) sin(7r$W/$o: 



dA(a 



Ml 



(i) 

Therefore, by setting $e = tvk, the qubit-resonator in- 
teraction can be turned off. When the interaction is on, 
$e can be tuned to compensate the difference of the fab- 
rication parameters and realize a homogeneous coupling 



(11) 
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g^- 1 = g. Then if each qubit is biased at the degeneracy in the interaction picture. Here, 9ij(n) = gWg^'T n /tJ — 
point = (n+ l/2)$ , the total Hamiltonian becomes g^g^2nn/uj 2 . Thus, at these times, the time evolu- 
tion is equivalent to that of a system of coupled qubits 

(i) (1) 

lc j with an interaction Hamiltonian of the form oc ct± <7i . 

Therefore, by choosing appropriate coupling pulse se- 
quences, an effective XX-coupling can be realized for mul- 
tiple qubits. This coupling can be utilized to construct a 
CNOT gate for two qubits 2 ^. If the couplings are homo- 
geneous for all qubits, i.e., g( % ' = g (for i = 1, -,N), 



(12) 

where fiW($£*^) = AM($i ) is the single qubit energy 
splitting. Comparing Eqs. (H]) and (fT!?)) , it is evident that 
the two Hamiltonians have the same structure: the in- 
teraction term commutes with the free term, and the 
interaction can be switched on and off. In the next sec- 
tion, we show how to generate a multi-qubit GHZ state 
by utilizing these features. 



II. GENERATION OF A GHZ STATE 

In the interaction picture, 

Hi(t) = J29 il) (^e luJt + ae- tuJt )a^ . (13) 

i 

bmce cri jS> , aax l) , aVi'' , 1} form a closed Lie Algebra, 
the time evolution operator in the interaction picture can 
be written in a factorized way as2£ 



tMo=n e-iA ii e 



(14) 



and Ui{t) satisfies 

i(^-U I (t))Uf 1 (t)=H I (t). (15) 

Solving this equation for the initial condition Aij(0) = 
B t (0) = D(0) = 0, we obtain 



Bi(t) = ^—{e- wt - 1) 



A ij (t) = 



o « o C0 / 1 

9 9 (J_^t_ 1) _ t 



U! \IU) 



oj \iu> 



(16) 
(17) 

(18) 



In the Schrodinger picture 

U s (t) = U (t)U I (t)=e- iua ' at l[, 



l Uj{t) . 



(19) 

Note that Bi (t) is a periodic function of time and van- 
ishes at t = T„ = 2irn/uj for integer n. At these instants 
of time, the time evolution operator takes the form 

U{T n ) = exp(-i^%(n)<7«4^)exp(-i£>(t)) , (20) 



An) = 9(n) = ^r2nn 



(21) 



Eq. (|20|) can be written as 

U(T n ) = exp(-iA9(n)J 2 )exp(i9(n)N)exp(-iD(t)) 



with J x =£i°-*V2. 

Suppose the initial state of the qubits is 



l*(0)) = 



i=l 



(22) 



(23) 



where \±) z denotes the eigenstates of a z , a z \±) z — 
±\±) z . This initial state can be prepared by biasing the 
qubits far away from the degeneracy point, letting them 
relax to the ground state and then biasing them back 
adiabatically. Starting from the initial state, under the 
time evolution described by Eq. (|22p . the state evolves 
into a GHZ stat o 28 i 29 (up to a global phase factor) 

1 / A x \ 

l*(T„)> = j= ( ® |->« + e^ N ^' 2 |+>« J , 

(24) 

if 9(n) = (1 + 4to)7t/8, where m is an arbitrary integer. 
A comparison with Eq. (f2"Tj) shows that the integers n 
and m are related by 



OJ 



Ag 2 I6g 2 ' 



(25) 



which is possible only if the (experimentally controllable) 
parameter g 2 /uj 2 is chosen to be 



1 + 4m 
16n 



(26) 



Since it is difficult in practice to realize g comparable 
to w, we assume m = 0. Hence Eq. (T2"6"l) determines the 
value n m ; n (typically larger than 1) which corresponds to 
the minimum preparation time of the GHZ state 



27rn n 



nuj 



(27) 



The optimal = 1 could be realized if it were 

possible to achieve g = lo/A. The same GHZ state is 
periodically generated at later times, with preparation 
time T p = T min (l + 4m). 
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For both types of qubits, g is proportional to i/uj (since 
g is proportional to 77, see Eqs. ([2]) and f4|). If we as- 
sume g = £-\/u), we obtain T mm = tt/8£, 2 . Therefore, the 
preparation time does not depend on to. Furthermore, 
the preparation time Eq. (|2T[) does not increase with the 
number of qubits. 

If the qubits evolve under the time evolution described 
by Eq. QTZty with 0(n) = (3 + 4to)tt/8, another N-qubit 
GHZ state is realized, 

/a n \ 

|ff(T B )> = ( (g) |->w + c -*<™>/» (g) |+)W J . 

(28) 

In the following discussion, we focus on the GHZ state 
Eq. (|24p since it can be prepared in a shorter time. 

The treatment discussed up to now is valid if the qubit 
number N is even. For odd N, the single-qubit rotation 
U' = cxp(— 7r J x /2) is needed in addition to the time evo- 
lution Eq. (|22p . The GHZ state that can be realized for 
odd N has the form 

/ JV AT \ 

|*(T, 1 )) = - 7 =f0|-)i 4) +e^/ 2 0|+)«J . (29) 

To conclude: one can prepare an iV-qubit GHZ state by 
turning on the qubit-resonator interaction for a specified 
time. 

For this GHZ state to be useful for quantum informa- 
tion processing, the preparation time has to be shorter 
than the quantum coherence time of the whole system. 
In general, a short preparation time results from a strong 
qubit-qubit coupling. However, this conflicts with the 
weak-coupling condition assumed in many schemes in or- 
der to utilize virtual photon excitation or the rotating- 
wave approximation. Our preparation scheme for the 
GHZ state is based on real excitations of the quantum 
bus. No weak-coupling condition is required here. In 
principle, it can be applied to the 'ultra-strong' coupling 
regime that the coupling strength between the quantum 
bus (i.e. the TLR) and the qubits is comparable to the 
free system energy spacing. Hence it is possible to im- 
plement GHZ state preparation in a very short time. To 
get an idea of the time scale under realistic experimental 
conditions, we now estimate the preparation time using 
typical experimental parameters. 

Assuming the mutual inductance between qubit and 
resonator AfW = 20 pH, the self inductance L = 
100 pH, and the resonator frequency to = 1 GHz leads to 
77 ~ 1.76 x 10~ 3 for both types of qubits. For charge 
qubits, we assume Ef = 14 GHz, fiW = 10 GHz, 
and that the bias during the coupling period satisfies 

sin(7r$i 4 ' ) /<I > o) = 0.8. This leads to a coupling strength 
of g = 19.71 MHz. For flux qubits, we assume a qubit 
frequency n(0 = 10 GHz, Ef = 345 GHz, ajf 3 = 0.42, 
the bias satisfies sin(7r$e /$o) = 0.71, and at this bias, 
dA/da = 112 GHz. Both 2a£° and 2af cos(7r$i l) /*o) 



should be within the interval (0.6,0.85) so that the cir- 
cuits can always work as flux qubits both with and with- 
out bias. This leads to g w 144 MHa 2 ^. The coupling 
strength is much stronger for flux qubits than charge 
qubits because of the direct magnetic coupling to the 
phase degree of freedom. 

Therefore the interaction time to realize a GHZ state 
is T m [ n = 1 jits for charge qubits and T m ; n = 19 ns for 
flux qubits. The preparation time for flux qubits is much 
shorter than the coherence time of the TLR which can be 
several hundred microseconds. The typical single-qubit 
coherence time at the degeneracy point is several mi- 
croseconds. Hence in principle, the scheme is able to 
prepare GHZ states for several tens of qubits. If the 
coupling strength can be further increased to the 'ultra- 
strong' regime in experiment, the preparation of a multi- 
qubit GHZ state can be comparable to the time of a single 
qubit operation. 



III. PREPARATION ERRORS 

From the above calculation, it is clear that the essential 
point to prepare the GHZ state is to control the length 
of the dc pulse to manipulate the flux bias 3>i^. In the 
beginning, the external magnetic flux 4>e is set to nn, all 
the qubits and the transmission line resonator are relaxed 
to their respective ground states. Then the interaction 
between the qubits and the resonator is turned on by 

(i) 

biasing $e away from nir to some appropriate value for 
a time T m i n . Finally, the interaction is switched off by 
setting = n7r again, and the multi-qubit GHZ state 
is realized. Note that all the qubit biases are modified 
during the preparation by the same pulse, therefore all 
the qubits can share one control line for the magnetic 
flux. To accomplish this operation, two practical issues 
have to be considered. 

The first one is the precision of the control of the pulse 
length to keep the error acceptable. If the pulse length is 
not exactly r m ; n , the state realized is not a GHZ state and 
this error can be evaluated by calculating the fidelit y 30 ! 31 
F(t) = Tr [pGHzPq(t)), where p q (t) is the reduced density 
matrix of the qubits and /9ghz is the density matrix of 
the 7V-qubit GHZ state. In Fig. [21 the blue curves show 
the fidelity of state preparation, the regime with fidelity 
larger than 90% is marked by two green dotted lines. To 
realize a preparation with above 90% fidelity, the time 
control of the pulse should be precise to around 2.5 ns in 
the four qubits case, which is possible in experiment. 

The second problem is the influence of the non-ideal 
pulse shape. In the above calculation, we have assumed 
that a perfect square pulse can be applied so that g^' is a 
constant during the preparation. However, in experiment 
the dc pulse generated always has a finite rise and fall 
time. Since the coupling strength g^ depends on the bias 

(i) 

flux <f>e , the modulation of the magnetic flux results in 
a time-dependent coupling strength g^ — g^(t). If g^ 
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FIG. 2: (Color online) Time dependence of the fidelity of the 
prepared GHZ state for two different initial resonator states: 
the ground state (blue line) and the thermal state (red line) 
in the case of (a) two qubits, (b) four qubits. The black dots 
indicate the time when the resonator and qubits are effec- 
tively decoupled. The green lines limit the regime in which 
the fidelity is larger than 90%. The following parameter val- 
ues were used: qubit frequency = 10 GHz, resonator 
frequency u = 1 GHz, coupling strength g = 144 MHz. The 
time is given in units of T m - m . 



varies slowly with time (compared with e~ lut ), the above 
discussions still hold except that the decoupling time T 
at which the qubit-resonator coupling can be canceled is 
shifted to satisfy 



-V ) (r)-ff w (o) = o 



(30) 



A GHZ state is prepared if 

f-T 



dt'{e^ t '[g^(t')g^{0)+g^(0)g^(t')} 
2< ? «(*')<? W (t')} (31) 



for all i, j. This means a GHZ state can be realized by dc 
pulses of finite bandwidth without introducing additional 
errors. 

Another systematic error appears because the param- 
eter g 2 /u> 2 cannot be controlled with arbitrary accuracy, 
i.e., Eq. (|26p will be satisfied only approximately. In 
experiment, $i is tuned to get the desired value of g; 
whereas to is fixed by the geometry of the device. Suppose 
the experimental inaccuracy leads to a modified value for 
the coupling strength, g{\ + 5), where 5 quantifies the 
magnitude of error. Hence the prepared state deviates 
from the GHZ state. The fidelity of the prepared state 
depends on 5 as 



where C r N — N\/(r\(N — r)\) is the binomial coefficient. 
This expression is valid for even N. For odd N, the 
fidelity turns out to be given by Eq. (|32p with N — > iV+1. 
Figure [3] shows that the fidelity decreases as the error in 
the coupling coefficient increases. In the case of a 4-qubit 
GHZ state, a fidelity of 98% can be achieved if the error 
in g is within 3%. However, as the number of qubits 
increases, the fidelity drops more rapidly. Hence a more 
precise control of the flux bias is required to realize many- 
qubit GHZ states. 



IV. ERROR CAUSED BY DECOHERENCE 

An important advantage of our proposal is that the 
state preparation is independent of the initial state of 
the resonator. In general, it is not easy to prepare the 
system to be exactly in the ground state. For example, 
at typical dilution fridge temperatures, say 50 mK, there 
is a non-negligible probability (30%) for the first excited 
state of a 1 GHz resonator to be occupied. This problem 
is less severe for the qubits since their energy scale is 
much higher. Therefore a scheme which is insensitive to 
the initial state is desirable. 

Figure [2 shows the fidelity of the prepared GHZ state 
for two different initial states of the resonator (the ground 
state and the thermal state at 50 mK) . Although the time 
evolutions are different in general, the fidelities at the de- 
coupling time T n (indicated by black dots in the figures) 
are the same. This can be explained from Eq. (fH)) , at 
times T n , only the first term of Eq. (fTij) is kept, i.e., the 
qubits and resonator are decoupled. No matter what the 
initial state of the resonator is, at these times, the res- 
onator has evolved back to its initial state. This means 
the GHZ state preparation is not influenced by the initial 
state, or in other words, the preparation is insensitive to 
the decoherence that occurred before the interaction was 
switched on. 

But the decoherence during the operation certainly 
changes the final output state. In general, environmental 
fluctuations induce both dephasing and relaxation to the 
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(32) 



FIG. 3: (Color online) Dependence of the fidelity F on the 
error of the coupling coefficient 8. The curves correspond to 
N = 2 (top), 4, 6, and 8 (bottom). 
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system. Since the qubits are all biased at the degener- 
acy point, the strong dephasing effect due to 1/f noise is 
largely suppressed. Thus we can use a master equation 
which only includes relaxation as damping instead of the 
unitary operator Eq. (fT4]) to fully characterize the time 
evolution 



p{t) = -i[H, p(t)] + C QP {t) + C RP {t) 



(33) 



where p(t) is the density matrix of the system (qubits + 
resonator) in the interaction picture and £r represents 
the decoherence of the resonator 



Cr P = 2 (-Nth + l)(2a / oa t - a)ap - pata) 
+ 2 -Nth (2a''' p° — aa) p — paa)) . 



(34) 



Here, k is the resonator decay rate and Nth = 
(exp(w/fcBT) — 1) _1 the average number of photons in 
the resonator. Finally, Cq represents the decoherence of 
the qubits 
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£qP = o (25--PO+ - pa+a~ - ct+5-_p) 



(35) 



where 7 is the qubit decay rate and a± are written in the 
diagonal basis of a x . The quality factor of a TLR can be 
as high as than 10 6 . The qubit Ti-time at the degeneracy 
point is several ps at most in present experiment. To be 
on the safe side, we assume for the resonator Q — 2 x 10 3 , 
k = 0.5 MHz, and for the qubit T\ = 100 ns, i.e., the de- 
cay rate is 7 = 10 MHz. Here we neglect excitations of 
the qubit since its energy spacing is much larger than the 
thermal fluctuation. To investigate the influence of deco- 
herence, we compare the fidelity to prepare a GHZ state 
with/without decoherence. The result is shown in Fig. [4] 
where the difference of the two fidelities AF = F — F<j 
(where is the fidelity in the presence of decoherence) 
is plotted as a function of time. The red dots mark the 
difference at the GHZ preparation times T p . Obviously, 
the error due to decoherence increases with time. As we 
analyzed in the previous section, the preparation time 
is much shorter than the decoherence time. Therefore 
the error is still quite small at the minimum preparation 
time Tinm (indicated by the first dot): the error caused 
by decoherence is around 3.7% in the 4-qubit case. 



V. DISCUSSION AND CONCLUSION 

In the above discussion, for simplicity, we assumed that 
the qubits only interact with a single mode of the res- 
onator. However, since we did not invoke the rotating 
wave approximation in our calculation, higher modes of 
the TL R 20 i 32 will also contribute to the coupling. There- 
fore, the interaction Eq. (flU)) should include a sum over 
multiple modes whose frequencies are below a cut-off 
lo c . The cut-off is determined by a number of practi- 
cal issues, e.g., by the superconducting gap, or the fact 
that the resonator is not strictly one-dimensional 20 . The 




FIG. 4: (Color online) Time dependence of the error due to 
decoherence. AF is the difference of the fidelity of the pre- 
pared GHZ state with/without environment decoherence for 
(a) 2 qubits (b) 4 qubits. The red dots mark the times at 
which the GHZ state is prepared. The following parameter 
values were used: qubit frequency O'*' = 10 GHz, resonator 
frequency w = 1 GHz, coupling strength g = 144 MHz, qubit 
decay rate 7 = 10 MHz, and resonator decay rate k = 0.5 
MHz. The time is given in units of T m i n . 



time evolution including higher modes is of the same 
form as Eq. (fl4|) but includes a product over all the 
relevant modes. Neglecting the small nonlinear effect 
due to output coupling, the frequencies of all higher 
modes are multiples of the frequency of the fundamental 
mode, LUn = nuj and u c — n c oj, all the coupling coef- 
ficients between the qubits and different modes of the 
TLR, B ifi (t) = igf{e- iu>Rt - 1)/uj a are still zero for 

t = T n = 2nir/to. Here g% = g^(u>n). Hence the only 
correction to our scheme is including a sum over all the 
relevant modes in the definition of Aij(t) in Eq. (fTT|) . 



™c (i) (j) 
Al3 (t) = V^-^i- 



n=l 



(36) 



For low excitation modes whose wave lengths are still 
much larger than the qubit dimension, the homogeneous 
coupling assumption still approximately valid, i.e., g\ = 
Pn- For example, considering h c = 10 for a 10 cm trans- 
mission line, around the center there is a 0.32 mm-long re- 
gion where the magnetic field varies within 5 %. The dis- 
tance between the center of two qubits is roughly 10pm. 
This means up to around 30 qubits are coupled to the 
resonator approximately homogeneously. One can also 
tunc to further compensate the slight inhomogeneity. 
The correction to the time evolution Eq. can be sim- 
ply written as 8(n) = (2irn)g 2 (n c /2)/ijj 2 . Therefore, the 
effect of the higher excitation modes actually amounts to 
increasing the coupling coefficient g — > gy/n c /2, which 
helps to reduce the operation time. 
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The electric field of the higher modes has little effect 
on the flux qubits but will change the voltage bias of the 
charge qubits and couple to their a z component. How- 
ever, at the degeneracy point, where the free Hamiltonian 
is proportional to a x , these coupling terms are rapidly os- 
cillating and are expected to have a small effect on the 
system. However, the situation is less advantageous than 
for flux qubits, and higher modes should be suppressed 
by choosing high fundamental mode frequencies in the 
charge-qubit case. 

For a small number of qubits, an (lumped) LC circuit 
can also be used as a quantum bu a 33 > 34 to generate a 
GHZ state by following our scheme. In this case, only 
one single mode contributes. 

In conclusion, we have proposed a scheme to prepare 
an 7V-qubit GHZ state in a system of superconducting 
qubits coupled by a transmission line resonator. We 
have analyzed the preparation scheme for both charge 
qubits and flux qubits. With this method, a multi-qubit 
GHZ state can be prepared within the quantum coher- 
ence time. In the case of flux qubits that is especially fa- 
vorable, the preparation time is two orders of magnitude 
shorter than the qubit coherence time. The preparation 



time can be reduced further if the coupling strength is 
increased to the ultra-strong coupling regime, where the 
coupling strength is comparable to the free qubit Hamil- 
tonian. The preparation scheme is insensitive to the ini- 
tial state of the resonator and robust to operation errors 
and decoherence. The coupling can be switched by dc 
pulses of finite rise and fall times without introducing 
additional errors. In addition, the scheme described in 
this paper utilizes a linear coupling which is intrinsically 
error-free if proper dc control is achieved. Due to all these 
advantages, this proposal could be a promising candidate 
for GHZ state generation in systems of superconducting 
qubits. 
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